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We discuss the possibility of realizing Weyl semimetal phase in the magnetically doped topological 
band insulators. When the magnetic moments are ferromagnetically polarized, we show that there 
are three phases in the system upon the competition between topological mass and magnetic mass: 
topological band insulator phase, Weyl semimetal phase, and trivial phase. We explicitly derive the 
low energy theory of Weyl points from the general continuum Hamiltonian of topological insulators 
near the Dirac point, e.g. k • p theory near F point for Bi2Se3. Furthermore, we introduce the mi- 
croscopic tight-binding model on the diamond lattice to describe the magnetically doped topological 
insulator, and we found the Weyl semimetal phase. We also discuss the dimensional cross-over of 
the Weyl semimetal phase to the anomalous Hall effect. In closing, we discuss the experimental 
situation for the Weyl semimetal phase. 



The discovery of topological band insulatorJiH^ has 
been bringing attention to the topological properties of 
the band insulators. These phases are insulating in the 
bulk but support metallic surface states which is pro- 
tected as far as the bulk gap remains open and time- 
reversal symmtry is respected. T he to pological prop- 
erties of topological band insulatorJ^IIII! are guaranteed 
due to the well-defined bulk gap. However, it's demon- 
strated recently that gapless Weyl semimetal phase can 
also have topological protection to opening up a gap in 
bulk, edge-bulk correspondence, and gapless anomalous 
Hall respons^^tiSl. 

In this paper, we present that topological insulator ma- 
terials can be turned into the Weyl semimetal phase upon 
doping the bulk with magnetic materials that is assumed 
to be ferromangetically ordered. Then, the magnetiza- 
tion mass and topological band gap start to compete and 
result three phases including topological Weyl semimetal 
phase. The Weyl semimetal phase shows up when the 
magnetization mass is strong enough to close the topo- 
logical band gap. This is striking in that in the previous 
studies on the effect magnetization (where magnetiza- 
tion mass is usually taken as the small 'perturbation'), 
it's thought not desirable to close the band gap as it 
will destroy the topological properties of the topologi- 
cal band insulator and turn the material into trivial one. 
We hereby show that the strong magnetization mass in 
three-dimensional topological band insulator materials 
can generate a new topological phase; Weyl semimetal 
phase. The other two phases are topological band insu- 
lator phase (characterized by the axion angle 6 ^ n) and 
trivial insulator phase {0 = 0). We show that the gapless 
Weyl phase is relatively easy to realize in the topological 
band insulator in experiment (See also the multilayered 
structure with topological insulators® which uses the 
surface bands of topological insulators and recently pro- 
posed material HgCr2Sei^, compared to the complex 
compound structure A2lr2 07*- ^ * ^ * i^ ^ where the possibility 
of realizing Weyl semimetal phase in condensed matter 
system is proposed. Our continuum theory for the Weyl 



semimetal phase is directly applicable to the well-known 
topological insulator Bi2Se;^. Our work would be per- 
haps most relevant to the material TlBi(Si_a;Sea;)2 where 
the gap could be tuned to the criticality in experimenlJSl. 
This continuum theory for the general low energy theory 
of the topological insulators is supported for the standard 
tight-binding model on the diamond lattice. Our discus- 
sion is then extended to the dimensional cross-over of the 
Weyl semimetal phase to the two-dimensional anomalous 
Hall effect. 



I. THEORETICAL PERSPECTIVE 

The universal low-energy physics of the topologi- 
cal insulator materials can be captured in the massive 
(isotropic) Dirac system 
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H = Vt'Y^ CJak'^ + Af (k)T^ (1) 

a=l 

where t'^ and are Pauli matrices. We treat as the 
'real' spin which is relevant in dealing with time-reversal 
symmetry. The Hamiltonian Eq.([T]) incorporates with 
the two discrete symmetries: the time-reversal symme- 
try and the inversion symmetry. Usually, the Dirac mass 
M has the momentum dependence, i.e., M — Mq + 6k^ 
which can be tuned to be in the topological band insu- 
lator phase if Mob < and trivial phase if Mob > 0. 
When we are in the topological insulator phas^, there is 
a single Dirac cone on the surface 
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iJ, =±W^(T'-ki (2) 

1=1 

which can support a quantized anomalous Hall effectP 
upon introducing time-reversal symmetry breaking gap 
^ ma^ on the surface. This effect attracts enumer- 
ous attentiorP^l^ which is related to "axion electromag- 
netism"EIl. The effect can happen when the magnetiza- 
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FIG. 1. Proposed phase diagram in terms of (m, M) for 6 < 
in Hamiltonian Eq.Q: (A),(B) Topological band gap M is 
larger than magnetization mass m, thus we are in the insu- 
lating phase. We have a topological band insulator (A) for 
M > and trivial insulator (B) for M < 0. (C),(D) Magne- 
tization mass is stronger than the topological band gap, and 
we have a Weyl semimetal phase. Note that if m — >■ —m, 
then two Weyl points change the sign of the chirality. At the 
transitions jm| — \M\, two Weyl points meet each other and 
result topological or trivial insulators. 



tion gap ~ m is safely smaller than Mq so that we are in 
the topological band insulator phase. 

We further ask what would be the effect of the magneti- 
zation mass if it turns large enough (due to the bulk mag- 
netic impurity doping). We assume that the magnetic 
material doped in the topological band insulator forms 
a ferromagnetic order to supply the necessary mass. If 
this magnetization is large enough, then strikingly it can 
generate a gapless Weyl semimetal phase. To study this, 
we introduce the magnetization mass in the form of 



Jm ■ ar"^ 



(3) 



to the Hamiltonian Eq.Q. Due to the isotropy in our 
model, we can further simplify Jm ~ mz (redefining 
Jm — > m) . As a whole, we will deal with the Hamiltonian 



H = vt''Y^ o-^k" -I- Af(k)r^ 



ma T 



(4) 



which supports three phases: topological band insula- 
tor with time-reversal symmetry broken surface, Weyl 
semimetal phase, and trivial phase (without losing gen- 
erality, we choose Mq > and m > 0) . We will show that 
the transition between two insulating phases and gapless 
Weyl semimetal phase happens at Mq = m. 

Upon with this general consideration, we start to dis- 
cuss the details of the Hamiltonian with the bulk mag- 
netization mass Eq.Q. 

- Energy spectrum of the band theory: The general be- 
havior of the spectrum of the Hamiltonian Eq. Q can be 



easily read off from the following commutation relations 
[a^T^,T^] = and [r^,T^cr^] ^ 0. Due to the commuta- 
tion [a^r'^, T^] — 0, we will have the competition between 
the magnetization mass m and the topological mass AIq. 
Explicitly, we have the energy spectrum E(k.) for general 
m and M (k) 

E{k) = ±^v^{k^^ + ky^) + {m± v/Af2 + v^kz^ (5) 

we find that there are two mass gaps A± where A+ = 
m -I- ^M'^ +v'^kz^ and A_ = m - ^/M^ + v'^k^^. We 
immediately see that the gap A_ can close if m > Mq 
at the momentum vk^ = ^^m^ — M(k)2 which are the 
positions of Weyl points in the momentum space. We 
now discuss two cases where m < A/q and m > Afg (See 
Figl). 

- Topological band and trivial insulator phases: Here 
we discuss the case m < Mq. To study these phases, 
we use the adiabatic argument, i.e., we are in the same 
topological class as far as the gap remains open at k = 
(we assume that k = point is the only point which pos- 
sibly closes the band gap) and the inversion symmetry is 
protected. As is always positive for any to, our main 
focus goes to A_ . The Dirac mass A_ (k = 0) = to — Mq 
which closes at to = AIq, signaling quantum phase tran- 
sition to the Weyl semimetal phase. When to < Mq, 
we need to further specify the mass M{'k) — Mq + bk^ 
and have a topological band insulator phase due to the 
band inversion if Mob < and a trivial insulator phase 
if Mgb > 0. When we are in the topological band insu- 
lator phase (and also the chemical potential lies in the 
surface magnetization gap ^ m < Mq), we have Chern- 
Simon effective surface theory with external electromag- 
netic gauge 



2m An 



(6) 



which is the quantum anom alous Hall effect with the 
quantized hall coefRcientP^ (Ta;y = ±6^/2/1 . 

- Gapless Weyl semimetal phase: We now proceed to 
the case to > Mq where the mass A_(k = 0) < 0. For 
the low-energy physics, we ignore the quadratic depen- 
dence on Ajf(k) = Mq + 0(k2) to clarify the physics of 
this phase. Then, we have two Dirac nodes in 0(k) where 
vkz = vKc,± — ±vKc = zty/m'^ — Mq"^. Near the critical 
points kz ~ i-ft'c, the low-energy theory of the Hamilto- 
nian Eq.Q reduces into (with qz^± — kz — Kc^±) 



H(k) - e(Ke 



u±qz,± vk- 
vk+ -u±qz,± 



(7) 



with anisotropic Dirac speed along z direction, u± 



This Hamiltonian is nothing but the Weyl 



fermion with the winding number ±1 localized at Kc,± 
(two Weyl points carry the chirality c± = sign(u±) = =f1 
at the Dirac point^). 

Due to this non-zero winding number, the Wwl Hamil- 
tonian Eq.Q has an anomalous Hall effeclPMl xy- 
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plane. Because each layer of the whole bulk participates 
in anomalous Hall effect, we define the anomalous Hall 
response per layer. The simple computation shows that 
the first chern number Ck^ — ±1 for kz S {—Kc, Kc) with 
Kc = \Kc^±\. Then, this would support the anomalous 
Hall conductance (when the chemical potential is exactly 
at the Dirac point) 

This effect can be defined up to modular e^//i due to the 
ambiguitj^^ in \Kc\ which is up to modular tt (half of 
the size of the reciprocal vector along z). 

When two Weyl points are brought to the Brillouin 
zone center and pair-annhilat^ (m — > Mq), we have 
topological insulator phase (characterized by axion an- 
gle 6 — ±7r) or trivial insulator phase (9 = 0). The two 
axionic response of the resulting insulator phases can be 
anticipated from the inversion symmetry encoded in Eq. 
^ and Eq. Q because the angle 9-^—9 under the 
inversion. By noting the following two reasonings (1) 
the axion angle 9 = ±7r {9 = 0) for topological insula- 
tor phase (trivial insulator phase) at to — ?> and (2) the 
angl e is t ied to topological structure of the band struc- 
ture liSMI and thus cannot be changed smoothly as far as 
the gap remains open, we conclude that the full region of 
topological band insulator (trivial insulator) phase in the 
phase diagram is characterized by 6* = ±7r (0 = 0). If the 
inversion symmetry is broken down, then the axion angle 
9 is no longer protected and expected to be continuously 
varying in terms of the parameters of the system. 

It's also interesting to note that we can realize quan- 
tized layered anomalous Hall effects if the two Weyl 
points are dragged {K^. ~ tt) to the Brillouin zone boun- 
day^ and pair-annhilated. Then each layer in the bulk 
carries the Hall conductance 

cr^y = — (9) 

and each layer forms integer quantum Hall states 
without magnetic field. From this consideration, we 
could remarkably obtain various interesting known topo- 
logical phases, topological band insulator (gapped), 
Weyl semimetal phase (critical), and three-dimensional 
anomalous Hall (gapped) phase, within the single effec- 
tive theory of the topological band insulator. In fact, we 
can go further if we notice that the same Hamiltonian 
can describe the two-dimensional topological band insu- 
lator. When the sample width of the three-dimensional 
topological insulator is reduced, then the material turns 
into the two-dimensional topological insulatoJ^il (with 
some oscillatory behavior as the function width). Fur- 
thermore, it's known that the two-dimensional topologi- 
cal insulator materials can be turned into the anomalous 
Hall effect upon doping of the magnetic 

materialPHlSl 

Hence, it's plausible to connect the three-dimensional 
Weyl semimetal phase to the two-dimensional anomalous 
Hall effect in the sense of the dimesional cross-over^ ^. 



From the above general discussion, we raise an inter- 
esting question: what will be the fate of the Weyl phase 
under reducing the sample width? To answer the ques- 
tion carefully, we consider that the sample width is re- 
duced along z, i.e., varies from some macroscopically 
large value (three-dimensional bulk limit) to a few layers 
(two-dimensional limit) and there are two cases we'd like 
to study: A. magnetization is orthogonal to the plane of 
the thin film (to = toz), B. magnetization is in the plane 
of the sample (to = mx) . We assume that the magnetiza- 
tion strength to is fixed while the sample width reduces. 
In the three-dimensional limit, we know that the both 
cases generate the Weyl semimetal phase if the magneti- 
zation TO is large enough. In the following discussion, we 
will find two different thin-film limits depending on the 
direction of to. When the sample width Lz is reduced 
enough, then kz starts to be discretized and each band 
can be labelled by G Z i.e., the electrons form the 
standing wave ~ smijiz'Kz/ Lz) in z. As the lowest en- 
ergy band Uz — 1 is enough to describe the physics for 
the thin film limit, we can simply replace kz by tt/Lz in 
the effective theory Eq.Q. We begin with the case B 
(to = mx) and the theory for the thin film is 

i? = Wr^ (Ja-\^+ M(k)T^ -I- TOCT^T° (10) 

a=l,2 

which has two Dirac points splitted in the momentum 
space along kx- Specifically, two Dirac points are at 
k^^c — ± ^"^^-^ i^ and we have the gapless phase. If 
the magnetization mass is not strong enough, we are in 
the trivial insulator phase as expected. For the case A 
(to = TOz), we would find the theory 

H^vt' ■ k'' + A/(k)r^ + ma'T° (11) 

a=l,2 

which is gapped everywhere. In fact, this is the model for 
the anomalous Hall effect if |to| > \M\ (otherwise, it's 
trivial), and the model supports the quantized anoma- 
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lous Hall response a-^y = 7;-. So, we conclude that the 
Weyl semimetal phase for the case A (m = mx) turns 
into the anomalous Hall effect if |m| > |M(k)| (or trivial 
insulator) as the sample width reduces. Note that while 
the sample width reduces, the bulk gap M(k) oscillated^ 
as the function of the width Lz and thus we expect to 
have the oscillation between anomalous Hall effect and 
trivial insulator. 

- Other commuting masses and possible gapless phases: 
The realization of Weyl semimetal phase in this paper 
rely on the commutation relation between mass/kinetic 
terms in Hamiltonian Eq.Q. In general, if two Dirac 
masses Q and M satisfy [Q, M] = 0, then two masses 
tend to kill each other. Given the topological mass 
M ~ and the kinetic term K ^ r'^a'^ in Eq. Q, 
we can classify the mass terms Q that gives the possi- 
ble gapless phases from the topological band insulator. 
The conditions are [Q, M] = and [if^jQ] 7^ for at 
least one of i = 1,2,3. Then these conditions include 
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Q = (rh ■ a)T^ which generates 'two' Weyl points, and 
Q = T^Gi which generates a gapless cy Under around the 
origin, e.g., k^"^ + ky^ ~ k^^ along k^ for Q = r'-^a'^ . 
However, the second case would be easily gapped out in 
perturbation. So the magnetization mass is the only mass 
that is capable of inducing the Weyl semimetal phase in 
competion with the topological band gap. 

II. ON THE REAL MATERIAL 

We begin with the basic model for Bi^^Sbi^^; with the 
magnetic mass due to the Zeeman coupling ~ Jm • S 
where m is the polarized magnetic moment in the sys- 
tem. Even though m can polarize in any direction, we 
further assume m to be parallel with z. This can be pre- 
pared by applying the external magnetic field in z when 
the magnetic order sets into the sample. So, we assume 
'mx,y « iTiz = TO and further ignore rrix^y We begin 
with the standard tight-binding model on the diamond 
lattic^ for the topological band insulator 



Hti = t ^ c-Cj -I- iX 



<hj> 



so Yl 

<<i,j» 



(12) 



where the second term encodes the spin-orbit coupling 
which connects the next nearest neighbors. The Hamil- 
tonian has the inversion symmetry and time reversal sym- 



metry. The spectrum of the Hamlitonian Eq.(12| is very 
well understood and is topological band insulator (with 
the distortion). Now, we introduce the magnetization 
over the sample and assume that the doped magnetic 
material approximately induces the uniform on-site Zee- 
man term 



(13) 



We move to the momentum space for the Hamiltonian 
Hti + Hz and find 



i/ = ^rf,(k)r + TOr34 



(14) 



where we followed the notation of Liang Fu et. alW We 
can diagnolize it to obtain the band structure 



E{\<i)^±^^dl + dl + [m±^dl + dl + dl)^ (15) 

Note the remarkable resemblance of the above energy 
spectrum Eq.(15l to the band spectrum Eq.([5]) which we 
obtained from the continuum model. As is indepen- 
dent of the momentum k, we don't need to solve for all 
momentum k but we can simply look at the points where 
the gap for Hti Eq.(12) is small and where (i3(k), (i4(k) 
and d^{k) go as ~ fkr+ O(k^). There are eight such 
points (r point, three X points, and four L points) in the 
Brillouin zone and we will have Weyl points out of these 
eight points. For example, we find that there are two 



Weyl points near the F point if to > m,c 
at 



kc = (0, 0, ± 



f2 



) 



v/to2 - (4t)2 



(16) 



where \At\ is the band gap at the F point. We find that 
this simple tight-binding model agrees with the contin- 
uum calculation, and hence we've provided that the mi- 
croscopic lattice model for the Weyl phase. 

- experimental preparation: In the experiment, topo- 
logical band gap M tends to be the largest energy scale 
in the system (~ 0.3eV for Bi2Se3^'') and this would act 
as the drawback to realize the Weyl semimetal phases in 
that the magnetization mass to should win against the 
band gap M. Hence it's desirable to reduce the band 
gap by bringing the system to near the critical point 
(M — )■ 0) between topological and trivial insulators. Re- 
cently, the mass gap of the topological insulator mate- 
rial TlBi(Si_a:Sea;)2 appears to be under the control and 
could be tuned to the criticality in experimentfSl. Though 
it's difficult to predict the energy scale for the magnetic 
mass, some theoretical calculations ('^ 5 - 50meV ) and 
experiments (~ O.leV) seems to be positiv^^^EH ^ see 
this physics (at least in the thin film limits). Due to the 
disorders, it might not be easy to detect the transition 
itself. However, if the magnetic mass is large enough, 
Weyl phase should appear. If there are more than one 
Dirac point in the low energy, there could be more than 
two Weyl points (depending on the precise form of the 
'representation' of the magnetization mass at the Dirac 
point). 

In terms of observation, the best way to detect Weyl 
semimetal phase would be perhaps ARPES experiment 
to see the Weyl points and the strange surface states. 
Transport experiment would be also applicable to detect 
anomalous Hall effects. Note that a^y ^ kc ^ (w? — 
M2)^/2 -vviiere M is tunable by controlling the width of 
the quantum well. Via this, <Jxy would be tunable in 
the thin film limit. For the reasobale prediction with the 
magnetization gap ^ 50meV in thin film limit and the 
optimal preparation of the sample Mq — > 0, we predict 
a^y ~ O(lOi) f7-icm-i. 

Note: Upon completion of this work, we became aware 
of a paper by A. A. Burkov ei. a/.'^ where the realization 
of Weyl semimetal phase by magnetic perturbation in the 
topological insulator (at criticality between topological 
insulator and normal insulator phases) is also considered. 



ACKNOWLEDGMENTS 

G. Y. C thanks Eun Gook Moon, Meng Cheng, Roger 
S.K. Mong, and Su Yang Xu for helpful comments. 
G. Y. C. also thanks Joel E. Moore for support and en- 
couragement during this work and acknowledges support 
from NSF DMR-0804413 and a KITP Graduate Fellow- 
ship. 



5 



M. Zahid Hasan and J. E. Moore, Annual Review of 



Condensed Matter Physics 2, 55 (2011), |arXiv:1011.5462 
I [co nd-mat.str-el] . 

^~J7 E. Moore and L. Balents, [Phys. Rev. B 75, 121306 



(2007)1 |arXi v:cond-mat/0607314 
L. Fu, C. L. Kane ' ^ ' 



ma 



W. Witczak-Krempa and Y. B. Kim, ArXiv e-prints 
(2011), [arXl^:1105. 6108 [cond-mat. str-e l]| 
3 H. Zhang,T;.-X. Liu, X.-L. QTxrDairz: Fang, and S.-C. 

Zhang, Nature Physics 5, 438 (2009), 
* S.-Y. Xu, Y. Xia, L. A. Wray, S. Jia, F. Meier, J. H. Dil, 
J. Osterwalder, B. Slomski, A. Bansil, H. Lin, R. J. Cava, 
and M. Z. Hasan, Scien ce 332, 560'(20ir)||arXiv:1104.4633 
Rev. B 79, 195322 (2009)[ |arXiv:cond-] [cond-mat.mes-hall], 



and E. J. Mele, Physical Review Letters 
106803 (2007), arXiv:cond-mat/0607699 
Phys, 



R^ Roy 
mat/060753ll' 

X.-L. Qi, T. L. Hughes, and S.-C. Zhang, Phys. Rev. B 
78, 195424 (2008), arXiv:0802.3537 [cond-mat.mes-halll 
M. Essin, J. E. Moore, and D. Vanderbilt, Physical 



Review Letters 102, 146805 (2009)„arXiv:0810.2998 [cond- 



mat.mes- iialil 



X. Wan, A. M. Turner, A. Vishwanath, and S. Y. 
Savrasov, Phys. Rev. B 83, 205101 (2011) , arXiv;1007.00l6l 
I [cond-mat.str-el] 

^ K.-Y. Yang, Y.-M. Lu, and Y. Ran , ArXiv e-prints (2011), 
|arXiv: 1105.2353 [cond-mat.str-el]] 
^ A. A. Burkov and L. Balents, ArXiv e-prints (2011), 



arXiv:1105.5138 [cond-mat.mes-hall] 

G. Xu, H. Weng, Z. Wang, X. Dai, and Z. Fang, ArXiv 
e-prints (2011), arXiv:1106.3125 [cond-mat.mes-hall] 
A. A. Burkov, M. D. Hook, and L. Balents, ArXiv e-prints 
(2011), arXiv: 1 110. 1089 [cond-mat.mes-hall] 



R. Yu, W. Zhang, H.-J. Zhang, S . -C. Zhang, X. Dai, and 
Z. Fang, Scie nce 329, 61 (2010)| |arXiv:1002.0946 [cond-| 
mat.mes-hall] 

""^7 L. Chen, J.-H. Chu, J. G. Analytis, Z. K. Liu, 
K. Igarashi, H.-H. Kuo, X. L. Qi, S. K. Mo, R. G. Moore, 
D. H. Lu, M. Hashimoto, T. Sasagawa, S. C. Zhang, I. R. 
Fisher, Z. Hussain, and Z. X. Shen, [Science 329, 659| 

[~(2010) 

'-' X.-L. Qi, R. Li, J. Zang, and S.-C. Zhang, Science 323j 
1184 (2009), arXiv:0811.1303 [cond-mat.mes-halll, 



L Garate and M. Franz, Physical Review Letters 104, 



1468 02 (2010), arXiv:0911.0106 [cond-mat.mtrl-sci] 

T. L. Hughes, E. Prodan, and B. A. Bernevig, ArXiv e- 

prints (2010), arXiv: 1010.4508 [cond-mat.mes-hall] 

^° A. M. Turner, Y. Zhang, and A. Vishwanath, Phys. Rev. 



I B 82, 24 1102 (2010) 

C.-X. Liu, H. Zhang, B. Yan, X.-L. Qi, T. Frauenheim, 



X. Dai, Z. Fang, and S.-C. Zhan g,|Phys. Rev. B 81, 041307| 
(2010)1 [arXiv:0908.3654 [cond-inat.mes-hall]| 



